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O ; Abstract 

Jh ■ Strongly magnetized nuclear matter within the context of Quantum Hadrodynamics (Walecka- 

■^ ■ model) and its extensions is investigated in this article. The magnetic field is coupled to the charge 

^f^ . and dipole moment of the baryons by including the appropriate terms in the Lagrangian density. 

The saturation density of magnetized, symmetric nuclear matter po{B) was calculated for magnetic 
fields of the order of 10^^ gauss. For the calculated range of po{B) the binding energy, symmetry 
^ . energy coefficient 04 and compressibility K of nuclear matter were also calculated. It is found 



that with an increasing magnetic field po{B) also increases, while the system becomes less bound. 
Furthermore, the depopulation of proton Landau levels leaves a distinct oscillatory imprint on K 



jy-s ■ and 04. The calculations were also performed for increased values of the baryon magnetic dipole 

f~^ , moment. By increasing the dipole moment strength po{B) is found to decrease, but the system 



becomes more tightly bound while the oscillatory behavior of K and 04 persists. 

PACS numbers: 26.60.Kp, 21.30.Fe, 21.60.Jz, 21.65. Cd 

Keywords: magnetized nuclear matter, relativistic mean-field, symmetry energy, compressibility 



* |jpwd@sun.ac.za] 
^ |fgs@sun.ac.za 



I. INTRODUCTION 

Investigating the properties of highly magnetized nuclear matter might shed some light 
on the behaviour of certain types of stars. All stars are magnetized, but the magnetic fields 
are typically not strong enough to influence the fundamental nuclear interactions in the 
stellar interior, except in the case of neutron stars. These stars are observed as rapidly 
rotating, stronglj^ radio-emitting objects called pulsars with magnetic fields of between 10® 
and 9 x 10^'^ G pL]. However, X-ray or 7-ray emitting pulsars are also observed. They are 
assumed to be highly magnetized pulsars with magnetic fields between 10^^ and 10^^ G and 
are called magnetars. For a review of their properties see Ref. |[2|]. 

Neutron stars are formed from the remnant core of a massive star that exploded in a 
core-collapse supernova. They are very dense, compact objects having the mass of 1.5 times 
that of the sun, but a radius of the order of 10 km. Therefore these stars have an average 
density similar to that of saturated nuclear matter. As a result it is generally assumed that 
these stars should consist, at least in part, out of cold dense nuclear matter. 

The current model of magnetars assumes that the star's magnetic field is formed in the 
interior of the progenitor through dynamo action |3|. Since the observed magnetic field 
strengths are those on the surface of the star, the interior field may be several orders of 
magnitude larger. As reported in Ref. [J], according to the scalar viral theorem these fields 
can be as large as 10^*^ G. Thus nuclear matter magnetized by a very strong, external mag- 
netic field could be expected in or near the crust of a magnetar. 

In this paper we investigate the impact that a very strong, external magnetic field has 
on the properties of cold, saturated symmetric nuclear matter for a range of the baryon 
magnetic dipole moment strengths. In particular we calculate the symmetric nuclear matter 
saturation density as a function of the magnetic field strength. This is done in order to 
establish the saturated configuration of nuclear matter for an increasing magnetic field: the 
conditions most likely found in the neutron star interior. Other properties that will be 
investigated are the binding energy, the compressibility as well as the symmetry energy, 
all at saturation. These characteristic properties have been established as indicators of the 



behavior of dense nuclear matter and are generally used to constrain nuclear matter models, 
such as Quantum Hadrodynamics (or QHD, which is also known as the Walecka- model) 5| 
and its extensions |6|. 

In QHD the nuclear interaction between the baryons stems from various mesons exchanges, 
where each exchange describes a different feature of the nuclear interaction. The free param- 

ike 



eters are the meson-baryon, as well as meson self-coupling, coupling constants. We will ma 
use of the QHDl [51, NL3 \^, FSUGold (or FSU) [8], as well as the more recent lU-FSU 



9| 



parameter sets. QHDl is one of the earliest parameter sets which include only the scalar 
sigma and vector omega mesons. NL3 also includes the isovector rho meson and introduces 
a self-coupling in the scalar meson field to improve the description of nuclear matter. FSU 
was parameterized to investigate, amongst others, the nuclei away from nuclear symmetry. 
It introduces a self-coupling in the vector field, as well as a coupling between the vector and 
isovector mesons. FSU and lU-FSU share the same couplings, but the latter was constrained 

n 

to also satisfy astrophysical requirements [9[. 

Despite the fact that all QHD parameter sets are constrained to reproduce the same nu- 
clear matter properties at nuclear saturation, they have very different behavior at densities 
above saturation |6|. In order to perform calculations the system is approximated using the 
relativistic mean-field (RMF) approximation. In the RMF approximation the meson fields 
operators are replaced by their ground state expectation values and become classical fields 

Various aspects of magnetized nuclear matter have already been investigated with QHD 
in the RMF approximation, most recently by Dong et al. [4']. In the latter work the density- 
dependence of the symmetry energy of magnetized matter was investigated with the FSU 
parameter set at various densities as well as proton and neutron ratios, while also adjust- 
ing some of the coupling strengths in the parameter set. The authors concluded that the 
parabolic isospin dependence on the energy per nucleon remains valid for strong magnetic 



fields. An overview of previous studies is also provided in Ref. j^ and references there-in. 



Q 



Casali et al. 10|] investigated the impact of magnetic fields of 10^^ and 10^^ G on the 



symmetry energy coefficient at densities below nuclear saturation using the NL3 and FSU 
QHD parameter sets. They found that the Landau levels gives rise to discontinuities in 
the symmetry energy and influence the composition of the neutron star crust. They also 
investigated the effect of including the coupling between the magnetic field and the baryon 
dipole moments and concluded that it will only be appreciated in very strong magnetic 
fields at subsaturation nuclear densities. However, it would appear that the contribution of 
the magnetic field was not included in the energy density of the neutron star matter. As 



pointed out by Broderick et al. in Ref. [ll|, the magnetic contribution should be included 



since it influences the equation of state and thus the composition of the matter. 



Broderick et al. llj were also the first to point out the importance and impact of in- 



cluding the coupling between the dipole moment of the baryons and the magnetic field, 
in addition to the coupling of the proton's charge to the magnetic field, in the description 
of magnetized matter. In Ref. Uj this coupling is referred to as the anomalous magnetic 
moment or "AMM" coupling with the coupling strength of the applicable baryon's magnetic 
dipole moment. We believe this to be a somewhat misleading term since baryons are not 
point particles like electrons, but have an internal structure of quarks and gluons. The 
anomalous contribution to the electron dipole moment arises from higher-order photon cou- 
plings to the electron charge. On the other hand, the baryon's dipole moment (partially or 
fully, depending on whether the baryon is charged or not) arises from the photons coupling 
to the baryon's charged internal structure. Thus these "anomalous" contributions to the 
baryon dipole moment are not higher order contributions of the electromagnetic coupling, 
but an expression of the fact that the baryons have internal structure. However, not to 
confuse the reader we will also adopt this naming convention. 

Due to this particular origin of the baryon magnetic dipole moment, more than anec- 
dotal evidence would suggest that it should not be constant under all conditions. One 
would expect that the internal baryon structure would be affected by the baryon density, 
especially at high densities. Since this structure is the origin of the dipole moment, by 
extension the baryon magnetic dipole moment would also be influenced. Experimental 
investigations have shown that the dipole moment of copper, which has one proton outside 
the closed Z=28 proton shell, with an even number of valence neutrons increase by about 



50% over a mass number range of 10 12|. Furthermore, Ryu et al. in Ref. 13| investi- 
gated the neutron star equation of state with density-dependent dipole moments for the 
baryon octet using the quark-meson couphng (QMC) models and extensions thereof. They 
found that that the neutron star equation of state is dependent on both the strength of 
the magnetic field as well as that of the baryon magnetic dipole moment. Unfortunately 
the density-dependence of the baryon dipole moment at high densities is not known, either 
experimentally or theoretically, which complicate this calculations. 



Since the purpose of this paper is to investigate the dependence of the saturation den- 
sity of symmetric nuclear matter po{B) on a very strong magnetic field, the following 
calculations were performed. Firstly, Po{B) was calculated as the density that minimizes 
the binding energy per nucleon for a fixed value of B. Then the compression modulus as 
well as the symmetry energy coefficient was calculated at these densities. Since the density 
dependence of the baryon magnetic dipole moment is not known, these calculations were 
repeated for a range of values of the magnetic dipole moment. The results and impact of 
the magnetic field as well as the adjustment of the dipole moment on dense matter systems 
such as neutron stars are discussed at the end of the paper. To start it all off, an overview 
of the formalism used is presented in the next section. 



II. FORMALISM 

The interacting part of the QHD RMF Lagrangian for magnetized matter is % |9| 
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is the isodoublet baryon field operator where p and n indicate the proton 



and neutron components, while (^q, Vq, and 60 indicate the scalar, vector, and isovector 
mesons. Furthermore, T3 is the isospin operator, A'^ = (0,0,Bx,0) with B = \B\, while 
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f [7'^,7''] are the generators of the Lorentz group [14]. The strength of the coupling 



between the baryons and A'^ is the baryon charge, qb 
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is the strength of the couphng between the magnetic and baryon fields (with units of the 
baryon magnetic dipole moment). Under normal conditions the proton dipole moment is 
2.793 /iAT while the neutron's is — 1.913 /ijv (expressed in units of the nuclear magneton fj,]y) 
15|. In Ref. [iG*] it is shown that the values of gn and gp equal to 
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reproduce the normal values of the dipole moments in the non-relativistic limit. To adjust 
the baryon dipole moments by a factor of x, gn and g^ should become [l6| 
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The second and ninth terms of Cmt dlD can be expanded to 



(3a) 
(3b) 
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where S 
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are the baryon spin-operators. The first term in (jl]) leads to the 



well-known Landau problem, where the single particle proton energy spectrum is quantized 
and consecutive levels differ in energy by factor of \qpB\. The last two terms represent the 



"AMM" coupling between 
baryons. As shown in Ref. 



;he magnetic field and the dipole density •0S'0 of each of the 



16J and references there-in, the magnetized proton spectrum is 

(5) 



e{k,, X,n) = \lkl+ Um*^ + 2\qpB\n + A^pfi j + g.V' + 
while the magnetized neutron spectrum is 

e(fc,A) = 
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ki + m*^ + XgnB + g^V - 



9pbo 



(6) 



where k\ = k^ + ky and m* = {m — gs4'o) is the reduced mass. When the magnetic field is 
orientated in the 2;-direction (in the 3-dimensional Landau problem) only the kz momentum 
contributes to the proton energy. With this particular choice of A^^ the proton spectrum is 
degenerate in ky while the particles are localized in the x-direction [16]. For neutrons the 
usual spherical symmetry of the ground state is broken by the magnetic field and replaced 
with cylindrical symmetry. 



In ([5]) and ([6]) A = ±1 and indicates the possible orientations of the baryon dipole mo- 
ment along z. Since the Pauli-matrices do not commute, the Hamiltonian of magnetized 
matter does not commute with the baryon spin operator S. Hence spin is not a good 
quantum number and we rather refer to the orientation of the dipole moment to distinguish 
between the two types of baryons which are usually referred to as spin up or down. There- 
fore, the "AMM" coupling induces a relative shift in the energy spectrum of both protons 
and neutrons with opposite values of A. 

For a system consisting of symmetric nuclear matter, like the one considered here, the 
density of protons and neutrons must be equal. Furthermore it is assumed to be at zero 
temperature and so the ground state will consist out of completely filled energy levels with 
energies up to the Fermi energy. The deformation of the neutron Fermi surface can be 
accounted for in the calculation of any neutron density, or ground state expectation value, 



by considering the directional momentum dependence of the Fermi energy 16 1 



Due to the Landau quantization, any magnetized proton density is expressed as a sum 
over the occupied Landau levels. The degeneracy of each proton level is incorporated by 
adding a pre-factor ^^j^-^- to the contribution of the level. The prefactor follows from a com- 
parison of the fundamental magnetic flux per particle to the total magnetic flux through 



the level (see Ref. 16|] and references there- in). 



Consequently the energy density of magnetized symmetric nuclear matter, including the 



contribution of the magnetic field, is Ref. 16 1 
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where is the Heaviside step function. A simplified expression for e can be found in Ref. 
1^, as well as references there- in. 
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FIG. 1. (Color online) Plots of the saturation density (left) normalized with regards to the B = 
values and the binding energy at saturation (right) for different QHD parameter sets with gi, = g\, ■ 



III. RESULTS 



Since the nuclear interaction is short ranged it can be saturated at very high densities. 
The saturation density is defined as the density at which the binding energy per nucleon 
is at a minimum (we follow the convention of the binding energy being negative) and the 
system is most strongly bound. The binding energy per nucleon is |5| 

e 
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— m 
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and most QHD parameter sets reproduce the saturation density at 0.148 fm ^ with a bind- 



ing energy of about —16.3 MeV 



9|. We calculated the saturation density as function of a 



,(0) 



strong magnetic field with Qb = g\, ■ The results are shown in Fig. [H 



From these plots we note that Eh{B) behaves very similarly for all parameter sets: as 
the magnetic field increases the system becomes less strongly bound. At magnetic field 
strengths of about 3 x 10^^ G, which corresponds to a density of between 1.2 and 1.4 times 
Po (depending on the parameter set), the system becomes unbound. Thus, albeit more 
weakly bound, the system can accommodate much denser matter as it becomes^ magnetized. 
The increase in po{B) with B was first noted by Chakrabarty et al. in Ref. [171]. However, 
in this paper the effect of the "AMM" coupling was not studied and the authors found that 
the system becomes more tightly bound with increasing B. 



While the various Eh{B) curves are similar, the rate by which po{B) increases differs. 



fiblMeV] 




,,Vo[MeV] 



450 










^^' 


400 






^ 






350 




^-"^ , 


-"' 


-' 


.'-•'' 




— -^^ 


1 






2 3 4 



B 

[lO'^G] "is 

m*[MeV] 
550 - ~ ~ 



500 




2.0 



2.5 



3.0 



B 



3.5 4.0 [10"G] 



450 




B 

[10"G] 



FIG. 2. (Color online) Nuclear matter observables at saturation as a function of B. Clockwise 
from top left we have the Fermi energy of the protons and neutrons (although only the neutron 
Fermi energy is shown since it does not noticeably differ from the proton values). Then next to 
it, n{X) labels the Landau level occupation of the protons for each value of A. Since gp = % < 
the difference between n(l) and n{—X) is 1. Furthermore, since B > the A = 1 protons have the 
lowest energies and thus occupies the most levels. In the bottom row m* and gvVo are plotted. 

As shown in Fig. [21 the various baryon Fermi energies fif, as well as the Landau occupation 
numbers n{X) do not differ significantly. Hence the Eh{B) curves would also be similar since, 
from Eq. ([8]) we deduce that the binding energy goes like the Fermi energy. However, m* 
and gvVo do vary significantly between the parameter sets. Since the underlying mechanisms 
are the same, these differences stem from the parameterization and must be the source of 
the variation in po{B). 



For QHDl gs(j)Q (from m* = m — gs4>o) and q^Vq have the largest values and so, since 
the Fermi energies are essentially the same, the QHDl densities will be the lowest with the 
fewest number of occupied Landau levels, as we see from the n(A)-plot. Applying the same 
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FIG. 3. (Color online) Compression modulus of magnetized nuclear matter for gb = db ^°'- 
various QHD parameter sets. On the right we show the same plot, but the values are normalized 
with regard to their values at po{B = 0). 

logic to the rest of the plots, it is deduced that the FSU parameter set will have the highest 
values of po{B). This is not surprising since in Ref. |9| (and references therein) it has been 
well established that FSU has the softest equation of state (EoS). Since the EoS relates 
the energy density to the pressure a stiffer EoS exhibits a more rapid increase in pressure 
with density. Thus it can accommodate higher (energy) densities with the (comparatively) 
smallest increase in the pressure of the system. 



Related to the stiffness of the EoS is the compression modulus K of nuclear matter which 
gives an indication of the compressibility of the matter. For magnetized matter the com 
pression modulus K{B) of the various QHD parameter sets were calculated using 



KB) 
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dp^ \p 



(9) 



P=Po{B) 



We note that Eq. ([9]) is essentially the derivative of the pressure as a function of density and 
thus an indication of the stiffness of the EoS. Since Po{B) increases with B, so too should 
K{B) since the higher the density the more incompressible the system becomes. From the 
plots of non-normalized values of K{B), on the left in Fig. |3l it is clear that QHDl has the 
stiffest EoS while FSU has the softest. On the right the normalized values of K{B), with 
respect to K{0), are shown. 

However, the increase in K{B) does not happen smoothly, but rather it fluctuates with 
an increasing amplitude (for small B these fluctuations are not visible on the scale of the 
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FIG. 4. (Color online) Analysis of NL3 K/Kq plot from Fig. [Sj of which a section is shown on 
the left. On the right we have the Landau level occupation. The vertical guidelines in both plots 
denote the same values of the magnetic field in both plots. 

figure). These fluctuations imply that as B increases the system varies the degree to which 
it is compressible. This behavior is related to the number of Landau levels occupied by the 

system. 

!n Fig. m the relation between the fluctuations of the NL3 K{B) and the depopulation 
19| of the occupied Landau levels are shown. Variable behavior in a magnetized system 
of charged particles is not unexpected, since in a (two dimensional) Quantum Hall system 
one observes dramatic variation of the conductivity due to the population of the Landau 
levels 20|. Consequently the influence of the Landau levels on K should not be entirely 



surprising. However, in this case the variability stems from the fact that the magnetic field 
influences the proton densities in two ways: through the degeneracy factor ^-^^ as well as 
the spacing between Landau levels which, from Eq. ([5]), goes like gpB. 

Therefore, as the magnetic field increases the degeneracy factor also increases which makes 
the system more compressible, but at the same time the energy gap between the Landau 
levels becomes bigger which makes the system more incompressible. Since the system is 
at it saturation density (hence in its lowest energy configuration), the Landau levels are 
depopulated one by one. As the occupation of the level with energy closest to the Fermi 
energy decreases the system becomes more compressible, since the particles from this level 
are absorbed by the Landau levels at lower energies. 



11 



a4(B) [MeV] 




70 






QHDl 


60 






NL3 
FSU 


50 







lU-FSU 


40 






_---"- 


30 


-_-_ r- 




20 








[10"G] 

FIG. 5. (Color online) The non-normalized (left) and normalized (right) values of 04 of magnetized 
nuclear matter for different QHD parameter sets. 

For neutrons B ^ only induces a relative shift in the energy of particles with differ- 
ent orientations of their dipole moments. For B > the A = — 1 neutrons {gn > 0) becomes 
the dominent particle since the A = 1 neutrons will flip their dipole moments in order to 
attain a lower energy A = — 1 state. This contributes to the increase in K{B), but has no 
influence on its fluctuations. 



To establish the asymmetric tendencies of a magnetized system the symmetry energy coef- 
ficient 04 is calculated. It gives an indication of whether symmetric or asymmetric matter 
is preferred (l8J. This coefficient is calculated at po{B) from 

1 d^ . 



04 



2 dt^ pb 

{N-Zf 



with t 



Pn- Pp 



(10) 



t=0 \ pb 

Since a^ is the coefficient of the '■^^'7^ term in the semi-empirical mass formula [18| the 
larger 04 becomes, the more symmetric matter is preferred in order to keep the energy at a 
minimum. 



The expression for 04 can be simplified, but care should be taken since the baryon Fermi 
energies are also dependent on the magnetic field. A simplified expression for 04 of magne- 
tized matter is given in Ref . [l6| . The results of the calculation are displayed in Fig. |3 



Similar to the other nuclear matter properties, 04 increases with B which indicates that 
more symmetric nuclear matter system is favored than when B = Q. The increase in 04 is 
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FIG. 6. (Color online) Comparison of the NL3 04 graph (left) and the difference in the proton 
densities (right). The vertical lines in both plots denote the same values of the magnetic field. 

intermittently interrupted by a slight decrease, signaling a preference for more asymmetric 
matter. Again this fluctuating nature of 04 is related to the depopulation of Landau levels. 
Its origin is best illustrated by considering the density of the two types of protons (A = ±1) 
and the difference in the proton densities App(A) = Pp(l) — Pp(— 1) is shown in Fig. |6l 

For neutrons, Ap„(A) increases linearly with B since the difference in the energy of neutrons 
with different orientations of their dipole moment is 2gnB. For protons, App(A) is not 
directly proportional to B since the difference in energy of protons with opposite signs of A 
depends on both 2gpB as well as the number of filled Landau levels for a given A. When 



A/9p(A) is at a local minimum An(A) = 1 2l(|, while at a local maximum (just before a 



Landau level depopulates) An{X) = 2 and one choice of A proton levels are preferentially 
filled. In between these points in App(A), a^ decreases and more asymmetric matter is 
preferred since the degenerate proton Landau levels can be filled at a lower energy cost than 
the (non-degenerate) neutron energy levels. 

Of further interest is to investigate the properties of magnetized saturated nuclear matter 
for adjusted values of the baryon magnetic dipole moments, which are adjusted according to 
Eqns. (I3bl) and ( l3al) . The strengths of the dipole moments were changed symmetrically so 
that the proton and neutron dipole moments are increased by the same factor, i.e. gt = ^Ogl 
means that the strength of both the dipole moments increased by a factor of 10. We observe 
that as gh is adjusted, the response of the different QHD parameter sets are very similar. 
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FIG. 7. (Color online) Plotted clockwise from top left are the saturation density, binding energy, 
K{B), and a^^B) for various values gi,, as multiples of g'^ . 

Hence only the results for the NL3 parameter set (an arbitrary choice) will be plotted. 
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In Fig. [7] the values of the different properties for gb equal to g^ , lOg^ , 20gl , and 
30gl are shown. We do not claim that these values of gb are necessarily feasible or attain- 
able for the plotted range of densities and magnetic fields, but rather that it illustrates the 
full spectrum of the possible behavior of magnetized matter under extreme conditions. 



It is observed that as gb increases the system becomes less dense but more bound. The 
density decreases since the gap between different dipole orientations of the protons and the 
neutrons, as well as the separation of the proton Landau levels, increase. Hence lower ener- 
gies are attained at the cost of the number of particles that the system can accommodate 
per unit volume. For very large values of gb the relative shift between the different A proton 
and neutron energy levels is so large that only one choice of A (the one with the lowest 
energy levels) is populated. 
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FIG. 8. (Color online) K{B) (left) and 04(5) (right) as function of po{B) for 5;, = g^ 



(0) 



The bottom row in Fig. [7] shows a4{B) and K{B) as (?;, varies for the NL3 parameter 
set. For both the fluctuating behavior persists but their respective increases becomes slow 
as gb increases, which mimics the manner in which po{B) increases with gf,. For gi, = 30gl at 
a magnetic field of the order of 3 x 10^^ G the last A = — 1 proton Landau level depopulates. 
At this point both K and a^ increases since the system becomes much more incompressible 
because only A = 1 Landau levels are now filled, while the system tends towards more 
symmetric matter as it can only accommodate A = 1 protons. 



IV. DISCUSSION 

It seems fitting that the results of this study are discussed within the context where the 
investigated conditions are most likely to occur: the magnetar crust and/or interior. 



From the results presented we see that for B < 10^^ G the saturation properties of nu- 
clear matter is not influenced significantly. However, as the value of gbB comes (relatively) 
close to that of the baryon (reduced) mass the system is undeniably influenced by the 
magnetic field [22|. To establish the scope of this effect it is useful to plot K{B) and 04(5) 
as a function of po{B), shown in Fig. [HI 



Shown there are the properties of saturated nuclear matter for an increasing magnetic 
field. This can be seen as a cross section of the neutron star crust, assuming that the 
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magnetic field increases with the depth of the crust. It is obvious that any sudden change in 
the magnetic field will drastically alter the properties of the saturated nuclear matter, since 
the configuration of the Landau levels will change. These changes have a direct impact on 
the compressibility of the system as well as the preferred mix of protons and neutrons. 

Highly magnetized neutron stars (magnetars) are known for various types of flares/bursting 
activities (see Ref. |2| for a review), some of which that are accompanied by glitches (sudden 
increase in the rotation frequency of the star) [23]. It has been reported by Woods et al. in 



Ref. [24!] that the flare in magnetar SGR 1900+14 was accompanied by the reconfiguration 
of the stellar magnetic field. Hence the changes in the stellar magnetic field, as well as the 
accompanying changes in the conditions within neutron star crust, may contribute to the ob- 
served properties of magnetars. If the reconfiguration of the magnetic field implies a change 
in the magnetic field strength throughout the crust, the changes in the compressibility can 
induce a compression wave in the stellar interior since the compressibility depends on the 
magnetic field in a nonlinear, fluctuating fashion. Since Eq. ([9]) is very close to the formula 
for the speed of sound in a medium c^, any such wave will travel at a different speeds at 
different densities in the stellar interior. If the changes in the magnetic field happen over 
a short enough timescale, differential stresses will build up in the crust due to the varying 
nature of K and c^. These differential stresses could be the trigger of a magnetar glitch. 

Turning to a^i^B) it is clear that any change in B implies a change in the preferred compo- 
sition of nuclear matter with regards to the ratio of protons and neutrons. If B undergoes a 
sudden change, energy could be released from either baryons flipping their dipole moments, 
or from the baryons within the crust undergoing inverse or normal /3-decay. These processes 
could be the source of some of the sudden bursts of radiation observed from magnetars. 

V. CONCLUSION 

We have presented some evidence that the composition and properties of nuclear matter 
depend on the strength of the magnetic field. The influence of the magnetic field thus cannot 
be ignored when studying the magnetar interior and could be the source of some of the 
observed properties of magnetars. We have also shown that the influence of the magnetic 
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field persists even if tlie strengtli of the baryon magnetic dipole moments changes. 

However, since the internal dynamics of the neutron star crust and, in particular, the 
time-dependence of the nuclear processes are not known it is difficult to model dynamical 
behavior. If the origin and behavior of the stellar magnetic field is known, then we might 
have a chance of calculating the infiuence of any changes on the crust. In a future pub- 
lication we consider whether a ferromagnetic phase might be present in the neutron star 
interior as the source of the magnetar magnetic field. 
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